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The nucleon momentum distribution riA(fc) for A —2, 3, 4, 16, and 40 nuclei is systematically ana- 
lyzed in terms of wave functions resulting from advanced solutions of the non-relativistic Schrodinger 
equation, obtained within different many-body approaches and different realistic nucleon-nucleon 
(NN) interactions. Particular attention is paid to the separation of the momentum distributions 
into the mean-field and short-range correlations (SRC) contributions. It is shown that, although at 
high values of the momentum k different approaches lead to some quantitative differences, these do 
not hinder the general conclusion that the high momentum behavior (k > 1.5 — 2 fm" 1 ) of all nuclei 
considered is a universal one, in that, to a large extent, it is independent of the mass number A, the 
used many-body approach, and the details of NN interaction. In order to analyze and understand 
the frequently addressed question concerning the relationships between the nucleus, UA(k), and the 
deuteron, no(k), momentum distributions, the spin (S)-isospin (T) structure of few-nucleon systems 
and complex nuclei is analyzed in terms of realistic NN interactions and many-body approaches. 
To this end the number of NN pairs in a given (ST) state, viz (ST) = (10), (00), (01), and (11), and 
the contribution of these states to the nucleon momentum distributions, are calculated. It is shown 
that, apart from the (00) state which has very small effects, all other spin-isospin states contribute 
to the momentum distribution in a wide range of momenta. It is shown that the T=0 and T=l mo- 
mentum distributions at k > 1.5 — 2 fm -1 exhibit, independently of A, very similar behaviors, which 
represents another evidence of the universal character of SRC. The ratio nA(k) /no(k) is analyzed 
in detail stressing that in the SRC region it always increases with the momentum and the origin 
of such an increase is discussed and elucidated. The relationships between the one- and two-body 
momentum distributions, considered in a previous paper, are discussed and clarified, pointing out 
the relevant role played by the center-of-mass motion of a correlated pair in the (10) state. Eventu- 
ally, the values of the the probability of high momentum components in nuclei and the per nucleon 
probability 02 of deuteron-like configurations in nuclei are calculated and critically discussed. 

PACS numbers: 21.30.Fe, 21.60.-n, 24.10.Cn, 25.30.-c 



I. INTRODUCTION 



Recent experiments on two nucleon knockout reactions at high values of the four-momentum transfer on carbon 
using protons, A(p,ppN)X , [lland electrons, A(e, e'pN)X, Q, as well as experiments on inclusive quasi-elastic (q.e.) 
electron scattering A(e,e')X 3H5[, have provided robust evidence on the long- hunted ground-state nucleon-nucleon 
(NN) short-range correlations (SRC), demonstrating that in both types of processes the projectile had interacted 
with a nucleon belonging to a correlated NN pair [6|. In exclusive experiments, where the knowledge of both the 
three- momentum transfer q, and the momentum of a fast detected proton p, allows one to reconstruct the momentum 
ki = p q = — Pmiss that the struck proton had before the interaction (provided the final state interaction (FSI) 
could be disregarded), it has been found [lj that in the region 1.4 < |p m i SS | < 2.8 fm" 1 the removal of a proton was 
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almost always accompanied by the emission of a nucleon N (mostly a neutron) carrying a momentum roughly equal 
to fei. At the same time, in the q.e. inclusive experiment A(e, e')X, the ratio of the cross section off a nucleus A 
to the cross section off the deuteron or 3 He in the region of the Bjorken scaling variable 1.5 < Xgj < 2 (the region 
of XBj where q.e. scattering off a correlated NNpair is expected to occur), exhibits a constant behavior, indicating 
that, in agreement with theoretical predictions [3], the virtual photon interacted with a nucleon of a correlated NN 
pair. The exclusive experiment, moreover, provided evidence not only on SRC in general, but also, in particular, of 
the dominance of proton-neutron (pn ) deuteron-like tensor correlations occurring in states (ST) = (10), as predicted 
by several realistic calculations [§l-llil] ; the experimental data [l[ have also provided information on the center-of-mass 
(cm.) momentum distribution of the correlated NN pair, finding agreement with predictions made long ago [l2^ . 
Recent reviews on experiments providing information on SRC and their theoretical interpretations can be found in 
Refs. [IH and [l4|. A detailed picture of SRC is however still limited to the 12 C nucleus, so that extension to other 
nuclei is necessary in order to have a general quantitative picture of SRC through the periodic Table. In a systematic 
study of SRC, particular attention should be given to the experimental and theoretical investigations of: (i) the relative 
and cm. momentum dependencies of SRC; (ii) their spin-isospin structure; (iii) the relative role of two- nucleon (2N) 
and three-nucleon (3N) correlations. We will discuss 3N SRC in a separate paper; here we concentrate on 2N SRC. 
In configuration space these can be defined as those deviations from the independent motion of two nucleons, moving 
in a mean field, when they approach relative distances r±2 — |ri — r2 1 = r < 1.2 — 1.5 /to; according to theoretical 
calculations, in this region, due to the very nature of the NN interaction (in particular to its central short-range 
repulsion and the tensor attraction in (ST)=(10) state), the two-body mean-field density is strongly suppressed and 
2N correlated motion dominates. The details of 2N SRC depend upon the spin-isospin state of the correlated NN pair, 
as well as upon the region of the nucleus one is considering, i.e. upon the cm. motion of the pair R = (ri + r 2 )/2. To 
investigate these details, one has to take advantage of the high momentum components generated by SRC that lead 
to peculiar configurations of the nuclear wave function in momentum space |7j. As a matter of fact, if nucleons "1" 
and "2" become strongly correlated at short distances, the momentum configurations, in the nucleus cm. frame, is 
characterized by ki ~ —fei, JCa-t. — Y2i=3 fej — 0> an d not by the mean field configuration ^2f =2 fei — — fei> i.e. when 
the high momentum nucleon is balanced by the entire A — 1 nucleons, each of them carrying an average momentum 
of the order ~ k/(A— 1). Thus, if a correlated nucleon with momentum k\ acquires a momentum q from an external 
probe, and is removed from the nucleus and detected with momentum p = k\ + q, the partner nucleon should be 
emitted with momentum fc 2 — —fei = Q P = P m iss- Such a qualitative picture is however strictly valid only if the 
cm. momentum of the correlated pair was zero before nucleon removal and, moreover, if the two correlated nucleons 
leave the nucleus without interacting between themselves and with the nucleus (A — 2). These effects have to be 
carefully evaluated when attempting to extract the momentum distribution from experimental cross sections. Within 
a mean field many-body approach, the main effect of SRC is to deplete the occupancy of single particle shell model 
states and to make the occupation of levels above the Fermi sea different from zero; this leads to a decrease of the 
momentum distribution at values of |fc| roughly less than the Fermi momentum hp, and to an increase of it, by orders 
of magnitude, with respect to the mean field distribution [Isl ]. In this context, it has been pointed out that even a 
low resolution measurement of the one-body momentum distribution at |fe| > 2 — 3 /to -1 , where mean- field effects are 
negligible, may provide precious information on SRC [l6| . Though the most direct way to investigate SRC would be 
via experiments that detect a pair of back-to-back nucleons in the final state, also experiments which are sensitive to 
the one-body momentum distributions could be very useful. We have analyzed two-nucleon momentum distributions 
in two previous papers [ToL [TT| ; here we concentrate on the one-nucleon momentum distribution nA(k), with the 
aim of clarifying some points concerning, particularly, its SRC and spin-isospin structures. We will quantitatively 
clarify to what extent the high momentum part of ^(fc) can be associated to the deuteron momentum distribution 
no(k). To this end we will show that such an association, which is qualitatively correct, has been motivated either 
from the results of approximate many-body approaches [17H20l |. or from just assuming it as an input in pioneering 
Monte Carlo many-body calculations of n^fc) [2l|. Recently, the momentum distributions of few- nucleon systems 
and complex nuclei have been calculated within sophisticated many-body approaches [T5L fill l2"2T - [34j , using modern 
realistic interactions (35Tj39j . For this reason it seems to us appropriate to update the situation concerning the 
relationship between the momentum distributions of a nucleus A, where all 2N spin-isospin states may contribute, 
and the momentum distribution of the deuteron, where only the state (ST) = (10) is present. Our paper is organized 
as follows. Our formalism, based upon proper spin-isospin dependent one- and two-body density matrices which allows 
one to calculate the various spin-isospin components of the nucleon momentum distribution n^(k), is presented in 
Section |nj In Section [IIII we: (i) provide some general definitions of the one- and two-body momentum distributions; 
(ii) illustrate the way SRC influence the momentum distribution; (iii) critically analyze the way the probability of 
SRC can be defined; (iv) present a systematic comparison of the momentum distributions of A=2, 3, 4, 16 and 40, 
nuclei resulting from different many-body calculations and NN interactions. The values of the calculated number 
of pairs in different spin-isospin states in a nucleus A, and the momentum distributions in these states are given in 
Section IIVI The comparison between the momentum distributions of complex nuclei and the deuteron momentum 
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distributions is illustrated in SectiorfVl In this Section the result of calculation of the probability of 2N correlations 
in nuclei is also presented. Finally, the Summary and Conclusions are given in Section [VTl 



II. SPIN-ISOSPIN DEPENDENT DENSITY MATRICES, MOMENTUM DISTRIBUTIONS AND 

SHORT-RANGE CORRELATIONS 

A. Nuclear ground-state wave function and spin-isospin dependent density matrices 

In this paper we consider the nuclear wave function of a nucleus with Z protons and N neutron (Z+N=A), resulting 
from the non-relativistic Schrodinger equation containing two- and three-body interactions, viz 

In Eq. ([l} tojv is the nucleon mass, and / and {x^}^ stand, respectively, for the set of quantum numbers of the state 
/, and the set of A generalized coordinates x^ = {r^, s^, tj}, with Sj and denoting the nucleon spin and isospin, and 

A 

Vi the position coordinates measured from the cm. of the nucleus (^2 r i — 0). Once tpf({xi}A) is known, various 

i=i 

density matrices pertaining to the nuclear ground state V"/=o = V'jMj can be calculated. For ease of presentation 
we will consider in what follows complex nuclei with zero total momentum J — in the ground state, and will use 
the notation ipQ = tp A . In this paper we are going to investigate the number of pairs in various spin-isospin states 
and the spin-isospin dependent two-body and one-body densities and momentum distributions. This requires the 
knowledge of two-body and one-body spin-isospin dependent density matrices, which can be obtained by introducing 
the spin-isospin projector operators Pff^ 1 = (3 + Ti ■ tj)/4 and and P^ =0 = (1 — Tj • i"j)/4, (with the same form for 
the spin operators). A list of the density matrices that we need in our calculations is given here-below: 
1. the non diagonal spin-isospin dependent two-body density matrix, viz 

A 

p^(r 1 ,r r ,r' 1 ,r' 2 )= / ^*(xi, x 2 . . . , x A ) £ ^ T \r u r 2 ; r[, r^)^(x' 1) x' 2) • • • , x' A ) [] dx^, (2) 
where the non diagonal two-body spin-isospin dependent density matrix operator is 

A A 

piP (n , r 2 ; r 2 , r' 2 ) = P? P% - r^ft - r 2 )*(r< - r'^ - r 2 ) J] S(r k - v' k ) J] S S3n ^ 5^ ^ , (3) 

k=£{i,j] n=l 

2. the half-diagonal two-body spin-isospin dependent density matrix, viz 



fi 2 

"2^E^ + £ i) 2(M)+ E MUM 

i—1 i<j i<j<k 



r A 

p^f (n , r 2 ; ri) = / ^*( Xl , x 2 . . . , x A ) £ p<f > (n , r 2 ; ^(x'l , x' 2) • • • , x' A ) J] ^ ( 4 ) 

l<j i— 1 

where 

A A 

tifHrur 2 ;r[)=P§PT j S(r i -r 1 )S(r j -r 2 )S(^-r' 1 ) J] *(r* - r' k ) J] 5 S3n ^ ^ ^ , (5) 

k^{t,j} n=l 

3. the diagonal two-body spin-isospin dependent density matrix, viz 

r A 

p^ (n , r 2 ) = / ^* (x lf x 2 . . . , SU) J2 Pif ] (n > ^ (xi , x 2 , • • • , x^) J] dx fc dx' fc , (6) 

J i<j k=l 

where 

^ T) (n, r 2 ) = P§ - r 1 )5(? j - r 2 ) J] S(r k - F' fc )5 S3fc s , <^ t , . (7) 



fc=i 
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The following relations between the various density matrices and their normalization should be stressed: 



(8) 



J P N s f> (n , r 2 ) dn dr 2 = jv&n , (9) 



where NrgT) is the number of nucleon pairs in state (ST), so that 

E = - do) 

(ST) 

As for the spin-isospin independent density matrices, they are normalized in the usual way, namely 
J p A (ri,r 2 )dr 1 dr 2 = A(A - l)/2, J ^(n, r 2 ; <) dr 2 = [(A - l)/2)] ^(n, r^), and / p A (ri)dri = A. Note that 
since the two-body state has to be antisymmetric, the possible (ST) states are: (ST) = (10), (01), L — even and 
(ST) = (11), (00), L = odd, where L is the relative orbital momentum of the pair. 

B. The spin-isospin independent and spin-isospin dependent two- and one- nucleon momentum distributions 

Having defined the spin-isospin dependent density matrices, we can introduce the two-body spin-isospin dependent 
momentum distribution of a pair of nucleons in state (ST), namely 

n^(k 1 ,k 2 )= JL_ J ^dr^n'^'e^-^^Oe^-^-n)^^,^;^,^). (11) 

By summing Eq. (1111) over T and S, the spin-isospin averaged two-nucleon momentum distribution is obtained 1 

n A (k 1 ,k 2 ) = £ ng^flfei.fca) = f dn dn' dr 2 dr 2 ' e lk ^ r ^) e ik '< r '- r '>) p A (ri, n; rj, r 2 ), (12) 

(ST) " 



where 



p^(n,r 2 ;ri,r- 2 ) = £ ^ 2 (n,r 2 ;n,r 2 ). (13) 

(ST) 

The two-body momentum distributions obey the following normalization 

J n^f(ki,k 2 ) dfex dfc 2 = J pj^ 2 (n,n) dn dr 2 = 7V(| T) (14) 



and 



J n A (ki,k 2 ) dk 2 dk 2 = J , 



A(A — 1) 

n A (ki,k 2 )dk 2 dk 2 = j p A (r\,r 2 )dr 1 dr 2 = . (15) 



In this paper we are interested in the various spin-isospin components n,g T s(ki) of the one-body momentum 
distribution of nucleon N± 

n ^ {kl) ^J^fJ J dridr'^A-^pA^A), (16) 



with normalization 



J n^ 1 (fe 1 )dfe 1 = J p^(n)dn=i- (17) 



In case of non isoscalar nuclei interference between different spin-isospin states may occur. Such a contribution in case of the three 
nucleon systems results to be negligible and will be omitted in the presentation. 
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where p^ 1 = pa/ A. More specifically, we have to find the spin-isospin dependent momentum distribution of a 
nucleon N\ that belongs to all possible N\N 2 pairs with given value of S and T. To this end, we need the two-nucleon 
momentum distribution of all pairs which contains nucleon Ni . Since the isotopic spin, unlike the spin, which is mixed 
by the tensor force, is a conserved quantity, we will first consider the two-body momentum distribution corresponding 
to a fixed value of T, i.e. the spin-isospin two-body momentum distribution, Eq. ([lip , summed over the spin S = 0, 1; 
this quantity will be denoted by n* liV2 (fci,fci), and, according to Pauli principle, we have 



~ , T L 2 (ki, k 2 ) 



\oo) 



(k 1 ,k 2 



N1N2 



(18) 



and 



^"(fci.Afc) - [n^(k 1 ,k 2 )+n i ^(k 1 ,k 2 ) 



N1N2 



n±n 2 1 



(19) 



where each of the four quantities n|g L ^ 2 (fci, k 2 ) is defined by Eq. (fTTj) . 

By integrating the two-body momentum distribution in isospin state T, we find the one-body momentum distribution 
of a nucleon N\ belonging to a pair with isospin T 



(iVijV 2 



(ki) = ^A J n^{k u k 2 )dk 2 



1 1 



ifei.(ri-ri) 



/ N 1 ,N 2 / 1 \ 7 



dn dr[ (20) 



with normalization 



J n^ N2 \k x )dk x = ^jj ^ lJV2 (fci,fc 2 )dfci dk 2 = 1. 



(21) 



where is the number of pairs N X N 2 with isospin T. In Eq. (f2"U| p^' N2 (ri, r[; r 2 ) is the half-diagonal two-body 
density matrix, Eq. Q , summed over the spin; it is the central quantity necessary to calculate two-body momentum 
distributions. The analogs of Eqs. (fT8|) and (fl~9| for the T-dependent one-body momentum distribution readily follow, 
namely 



'(00) 



(22) 



4=f 2) (fc0 



N X N 2 
l (01) 



(23) 



In order to obtain an explicit equation for the momentum distribution we have to know the weights of a given 
isospin state in nucleus A. Let us consider the proton distribution, which gets contributions from pn and pp pairs; 
the former can be in T — and T = 1 states, whereas the latter can only be in T=l state. We have therefore to find 
the weight of T — and T = 1 pn pairs in nucleus A, since the weight of pp pairs in T = 1 state is one. The total 
number of pn pair in T=0 state in nucleus A with isospin Ta is (see Ref. |44| ) 



A^Lo = N& + N& = I [A(A + 2) - 4T A (T A + 1)] . 



(24) 



Dividing Eq. (|24l) by the total number of pn pairs, NZ, we find the weight Wj,_ of a pn pair in nucleus A, namely 



w; 



pn 
T=0 



8ZN 



[A(A + 2)-4T A (T A + l)], 



(25) 



with the weight of a pn pair in T = 1 given by wtp =1 = 1 — w^_ . Thus, we obtain the momentum distribution of 
nucleon iVi in terms of the explicit T = 0, 1 contributions 



1 



S± — X L 

which is correctly normalized to one since w^ =Q 



-{Z [< n =0 n^(fex) + v$U 4^1(^1)] +(Z- 1)4^^1)}, 



(26) 



u T=l 



1 and both n 



(N 1 N 2 ) 



k\) are normalized to one and where 



J n^P =0 {k\)dki — J n^P =l {ki) dk\ — J nf?_ 1 (fei) dk\ = 1. 



(27) 
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Note that using Eq. (|25l) in Eq. (|26p an even simpler equation is obtained for isoscalar nuclei, namely 



1 



A — l 



2 n^ (fci) + 3 ^ 2 n^ x (fei) 



Using Eqs. ^ and 



4 i=uv-y ■ 4 

we can write 

ru^i)^ = ni 10) (fei) +ni 00) (fei) + r# 1} (fci) +n^ 1) (fe 1 ), 



(28) 



(29) 



where all A-dependent coefficients are incorporated in the proper n^ T ^ . 
The calculation of the quantities n A T \k\) will be presented in Section HV 



III. NN INTERACTIONS, MANY-BODY APPROACHES, NUCLEON MOMENTUM DISTRIBUTIONS 

AND SRC 

We will now address the question concerning the content of SRC in the nuclear wave function, in particular the 
question concerning the definition of the probability of two-nucleon SRC, for, here, a certain degree of ambiguity may 
easily arise. The ground-state wave function V'jMj * s the solution of the many-body Schrodinger equation and it 
describes both mean- field and correlated motions. The latter includes both long- and short-range correlations; long- 
range correlations manifest themselves mostly in open shell nuclei, making partially occupied states which are empty 
in a simple independent particle model, with small effects on high momentum components; short-range correlations, 
on the opposite, generate high virtual particle-hole excitations even in closed-shell nuclei, and strongly affect the high 
momentum content of the wave function. Therefore, assuming that the momentum distributions could be extracted 
from some experimental data, we have to figure out a clear cut way to disentangle the momentum content generated 
by the mean-field from the one arising from SRC. To this end let us use the following procedure If we denote 
by {\ipf _1 >} the complete set of plane waves and eigenfunctions of the A—l Hamiltonian of the A — l nucleus, 
containing the same interaction as the Hamiltonian which generated the ground-state wave function and use the 
completeness relation 



E 



= i 



into Eq. (116l) . it is easy to see that the one-nucleon momentum distribution becomes [4C 

n^(k 1 )=n^(k 1 )+n^(k 1 ), 



where 



(2tt 



1 ]T ||e lfe -^Mr 1 | X | ffi 41- 1 )*(r 2 ...r A )V, A (r 1 ,r 2 ...r A )n^ i 



/=0,<71 



and 



— t 



f e ik ^dn f x\ a y f % 1) *{r 2 ...r A )^{r l ,r 2 ...r A )X[dr 

J J 2 i=2 



(30) 



(31) 



(32) 



(33) 



In the last equation the sum over "/" stands also for an integral over the continuum energy states, which are present 
in Eq. (|30[) . We see that the momentum distribution can be expressed through the overlap integrals between the 
ground-state wave function of nucleus of A and the wave function V'/ A_1 ' ) of the state / of nucleus (A — l). The 
squared modulus of the overlap integral represents the weight of the ground and excited states of (^4 — 1) in the ground 
state of A, so that the quantities 



and 



V Nl = 
1 gr 



ex 



'yr 



ki)dki 



n^(fci)d*!i. 



(34) 



(35) 
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with 



r£ + r£ = i, (36) 

can be associated to the lack of ground-state correlations (V^. 1 ) and to the presence of them (P^ 1 ). The separation 
of the momentum distributions in n^, 1 and n^. 1 is particularly useful in the case of A = 3, 4 systems, i.e. when the 
excited states of (A — 1) are in the continuum. In the case of a complex nucleus, where many discrete hole excited 
states are present, it is more convenient to use another representation where the particle-hole structure of the realistic 
solutions of Eq. ([l} is explicitly exhibited, namely 

ipa(rx,r 2 ...r A ) = c Q <5>o p0h (r 1 ,r 2 ...r A )+ c 2 <f> 2 A p2h (r 1 ,r 2 ...r A ) + ... (37) 

In Eq. (07} , oh is a Slater determinant describing the mean field motion of A nucleons occupying all state below 
the Fermi level, ^ 2p2 h describes 2p-2h excitations due to SRC, and the dots include higher order p-h excitations. The 
modulus squared of the various expansion coefficients Ci is nothing but the probability to have np — nh excitations 
in the ground state wave function. In particular \c 2 \ 2 = a 2 will determine the amount of ground-state SRC. Within 
such a representation, one can write [4l| 

n^(k 1 ) = n^(k 1 )+n^(k 1 ), (38) 



where 



o 1 (ki) = j^ E £ \je ik ^d ri j X \^- 1 ^(r 2 ...r A )^(rur 2 ...r A )fldr i \ 2 (39) 



f<F,<r 



and 



^0 = 7^3 E | / * ik *' ri drx I x\ a J?- 1) \r 2 ...r A )^{r l ,r 2 ...r A )Y[dr i ] \ (40) 



S>F,a 



where the summation over / in Eq. (|39|) includes all the discrete shell-model levels below the Fermi level (F) in 
(A — 1) ("hole states" of A), and in Eq. (|4T)]) it includes all the discrete and continuum states above the Fermi sea 
created by SRC. In a fully uncorrelated mean- field approach, we have 



n A h (k 1 ) = n^(k 1 )=J2\Mki)\ 2 ; n^(h) = 0, (41) 

a<F 

and the the analogs of Eqs. (jM)) and ([551) arc 

K 1 = J n^(kx)dk! (42) 

V? 1 = J n^ik^dkx (43) 

with 

■pN 1+ -pN 1 = 1 ( 44 ) 

The quantities V^ 1 and Vq 1 yield, respectively, the probability to find a mean-field and a correlated nucleon in the 
range < k\ < oo. Therefore they can be assumed as the mean-field and SRC total probabilities. It is clear that both 
low and high momentum components contributes to mean field and correlated momentum distributions, but we shall 
see, as expected, that n^ 1 (n^ 1 ) gets contribution mainly from low (high) momentum components. Assuming that 
Uq 1 and n 1 ^ 1 could experimentally be obtained, it might well be that only a limited range of momenta is available 
experimentally, in which case it is useful to define the partial probabilities 
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n O(l)( fc l) fc l dfa 



(45) 



i.e. the probability to observe a mean field (correlated) nucleon with momentum in the range k x <ki <k^. Although 
we do not discuss in this paper how the momentum distribution could in principle be extracted from the experimental 
data, we would like nevertheless to briefly comment on this point. It is clear from the very definition of the momentum 
distribution, that in order to obtain information on it one has to figure out an experiment in which a nucleon is struck 
from a nucleus A and the nucleus (A — 1) is left in a well defined energy state. To fully understand the point, it is 
useful to introduce the nucleon Spectral function, i.e. the following quantity 



S^{k x ,E) =< ^\al uax 8{E -H + E A )a kl>rTl \^ > 



E 



■r A )Y[d ri S(E-E 



J A-1 



Ea)) 



= S^(k 1 ,E ) + S 1 ^(k 1 ,E), 



(46) 
(47) 
(48) 



where in Eq. (|4*5|) . which has been obtained from Eq. (j4"Tf using the completeness relation (Eq. ([50])). aj^ (a^) * s a 

creation (annihilation) operator, E A = M A — M A -i — ttin , and E = E A + E A _ X is the nucleon removal energy. One 
has 



S^(k 1 ,E) = 



f<F 



e iki-r ldri 



, (A-l) * i 



r A )^(r 1 ,r 2 ...r A ) J]rfr, 6{E - E A _ X - E A ), (49) 



i=i 



A 

S^(k 1 ,E) = J2\ [e ik ^d ri f x{ ai 4 A ' 1) *(r 2 ...r A )^(r 1 ,r 2 ...r A )l[dr i \ 2 S(E-E{_ 1 -E A ). (50) 

The spectral function represents the probability that, after particle " 1" is adiabatically removed from the bound state 
and placed in the continuum, the nucleus (A — 1) remains in the state E^_ v The relation between the spectral 
function and the momentum distribution is given by the momentum sum rule 



J S A il {k u E)dE = n A h (k 1 



(51) 



The partial and full momentum distributions can therefore be obtained in principle by detecting the final nuclear 
system {A — 1) in correspondence of / < F and f > F. The exclusive processes A(e, e'N)(A — 1)/ in Plane Wave 
Impulse Approximation (PWIA) depends directly upon <S(ki, E). Thus by performing these types of experiments in a 
wide range of excitation energies of the final (A — I) nucleus and by performing the integration over E the momentum 
distributions can be obtained. FSI make the cross section to deviate from the PWIA, and, moreover, for a complex 
nucleus the sum over the entire continuum spectrum of (A — 1) is difficult, if not impossible, to perform. In case of 
few-body systems this difficulty can be overcome, since the number of possible final states is strongly reduced and, as 
a matter of fact, experimental information of n gr and n ex for 3 He and 4 He is already available [42 



43] 



We reiterate that the aim of this paper is the theoretical investigation of some general properties of momentum 
distributions, concerning in particular their SRC and spin-isospin structures. To this end for A=3 and 4 "exact" wave 
functions obtained, either by a direct solution of the Schrodinger equation or by variational procedures, will be used, 
whereas for complex nuclei momentum distributions obtained from various methods, ranging from the Brueckner- 
Bethe-Goldstone approach to the cluster expansion techniques will be adopted. In the next Section the momentum 
distributions of several nuclei will be presented and the values of the quantity T^Q^ikf) (Eq. (1431) ) will be given. 



A. The momentum distributions of few- nucleon systems and complex nuclei 

In this Section the momentum distributions of 3 H, 3 He, 4 He, 16 and 40 Ca, calculated within different approaches 
and using various two-nucleon interactions, will be presented. The full momentum distributions are shown in Figs. 
[T][B1 whereas their separation into the mean- field and correlation contributions, according to Eqs. ([52]) . (|53|) . (|39p and 
(|40|. are presented in Figs. [7lfT0l Note that from now on the notation k = ki and k = \ki\ will be used. 
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1. The momentum distributions of 2 H 



The momentum distributions of 2 H obtained by solving exactly the Schrodinger equation is crucial for our analysis. 
It is presented in Fig. [TJ where it can be seen that, apart from the RSC interaction [35j], the Paris interaction (3(| 
and the family of Argonne interactions AV8' 37 1, AV14 [38|, and AV18 39;!], provide essentially the same result. All 
these potential exhibit a strong short range repulsion which gives rise to a strong suppression of the deuteron wave 
function at internucleon separation r = |ri — r2 1 ^ 1-5 fm . This, together with the effects from the tensor force, 
generate high momentum components in the momentum distribution. 



2. The momentum distributions of 3 H and 3 He 

As already stated in Section IIII1 the three- and four-nucleon systems 3 H, 3 He and 4 He are very important, in 
that n gr and n ex have been explicitly calculated within accurate few-body techniques. Moreover, being 3 He and 3 H 
non-isoscalar nuclei, their proton and neutron distributions are different. As a matter of fact, in 3 He the proton 
momentum distribution is given by 

nl{k) = nl r {k) + nl x {k) (52) 

and the neutron distribution, due to the absence of a two-body bound state in the final state (cf. Eq. (|33|) ) , is given 
by 

nm = (53) 

In the above equations, n p gr {k) is the Fourier transform of the overlap between the ground-state wave functions of 
3 He and 2 H (cf. Eq. ([321 ). and nP x (k) the Fourier transform of the overlap between the ground-state wave function 
of 3 He and the continuum state of pn pair (cf. Eq. ([33])). Thanks to isospin invariance, Eqs. ([52|) and (|53|) represent, 
respectively, the proton and neutron momentum distributions in 3 H. The proton and neutron momentum distributions 
in 3 He resulting from Faddeev and variational calculations in correspondence of the AV18 interaction are shown in 
Fig. [2] It can indeed be seen that they are different, with the former strongly differing from the deuteron momentum 
distributions. The origin of such a difference will be discussed in detail in Section ITVl 



3. The momentum distributions of 4 He 



The nucleus He is the lightest isoscalar nucleus, with identical proton and neutron momentum distributions. These 
have been calculated in Ref. 25[ within the approach of Ref. 24| using the AV8' interaction. They are compared in 
Fig. [3] with the results of the variational Monte Carlo method performed with the AV14 interaction [22}. 



4- The momentum distributions of O and Ca 



The momentum distributions of complex nuclei is by far more complicated to calculate with the same accuracy 
attained in the case of three- and four-nucleon systems. Nonetheless, several calculations for 16 have been performed 
within different approaches and using various NN interactions, namely with the RSC potential |35|. in Ref. [17H20| . 
with the AV8' potential [37j], in Ref. [H] and Ref. [Hj], and with the AV14 potential [H] in Ref. [12] ; the various 
methods that have been used are the unitary operator approac h fl7j , th e Brueckner-Bethe-Goldstone approach [l8l[20|. 
the cluster expansion approach truncated at different orders |19l |3l|, the Fermion-Hypernetted-Chain method [28j . 
and the variational Monte Carlo correlated approach 22]. The various results are compared in Fig. |U As for 40 Ca 
two available results obtained with the V8' interactions are shown in Fig. [5] 



5. The A-dependence of the momentum distributions 

The momentum distributions of the considered nuclei obtained with the V8' interaction (the AV18 in the 2 H and 
3 He cases), are compared in Fig. [51 The general features that emerge from such a comparison can be summarized 
as follows: i) at low values of the momentum k = |ki| the shape of riAik) is determined by the asymptotic behavior 
of the wave function of the least bound nucleon, and therefore it is very different for different nuclei, ii) in the high 
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momentum region (k ^ 1.5 — 2 fm ) a qualitative similarity between the momentum distributions of deuteron and 
heavier nuclei can be observed. In what follows we will show that in this region riA(k) is dominated by the correlated 
part of the distributions, namely n ex and ni(k), and that the similarity between deuteron and complex nuclei is 
only a qualitative one, with the high momentum behavior of riA(k) being governed by the the various spin-isospin 
components contributing to n^fc), and not only by the deuteron-like state (ST) = (10). 



6. The mean field and SRC contributions to the momentum distributions 



The separation of the momentum distribution according to Eqs. (f5T|) and is shown in Figs. [71ITU1 It can be 
seen that: (i) in the region k < 1.5 — 2.0 fm~ x SRC reduce the mean field distribution without practically changing 
its shape, the effect being due to the decrease of the occupation probability of the shell model states below the Fermi 
level; (ii) in the region k > 2.0 fm~ l the momentum distribution are entirely exhausted by SRC. Having at disposal 
both n gr (k) (no(fc)) and n ex (k) (n\(k)) the probabilities given by Eqs. (|55l) . (|4*2"j) and (|4*5|) can be calculated. 

These are listed in Table HI whereas the partial probabilities defined by Eq. PS")) are listed in Table |TTJ 



B. Summary of Section [TT1 



From what exhibited in the present Section, some general features of the momentum distributions can be identified, 
which are, to a large extent, independent of the many-body approach and the two-nucleon interaction used in the 
calculations, namely: i) at k > 2 fm" 1 the momentum distributions of both few-nucleon and complex nuclei quali- 
tatively resemble the deuteron momentum distributions; (ii) in the region of high momenta, the realistic momentum 
distributions of complex nuclei overwhelms the mean-field distributions by several orders of magnitude; (iii) whereas 
for few-nucleon systems the method of calculations is very well established, for complex nuclei different methods and 
potentials provide at high momenta values of the distributions which can differ up to a factor of two, and it is not yet 
clear to which extent such a difference should be ascribed to the different potentials or to the different methods; (iv) 
the momentum distribution extracted from A(e, e'p)X and from the y-scaling analysis of inclusive A(e, e')X scattering 
[r| agree with many-body calculations; although the errors of the extracted momentum distributions are very large 
at high momenta, they are much smaller than the difference between correlated and mean-field distributions, with the 
latter being totally inadequate to predict high momentum components. In what follows our analysis of the momentum 
distributions will continue using the most advanced available calculation methods and two-nucleon interactions. In 
order to understand the microscopic origin of the correlated part of riA(k), we will analyze in the next Section its 
spin-isospin structure. 



IV. THE SPIN-ISOSPIN STRUCTURE OF THE NUCLEON MOMENTUM DISTRIBUTION AND SRC 



The spin-isospin structure of SRC is a fundamental quantity since it reflects the details of the NN interaction in 
the medium. It is therefore important to investigate how such a structure can affect various quantities which are 
related to SRC, like, e.g., the nucleon momentum distributions. In Ref. [23| a detailed analysis pertaining to few- 
nucleon systems has been presented of the various (ST) channels contributions to the relative momentum distribution 
n {ST) 2 (^l) ^2), Eq. (fTTj) , integrated over the cm. momentum, namely 

n (ST)fo re i) ~ J n (ST) (^i> ^2) <fK c TO = J grpj (k re i , K c rn ) dK. crn , , (54) 

whereas in Ref. [HI the dependence of the two-body momentum distribution n^^f 2 (fei, fe 2 ) = n^g^f) (k re i, K c , m ,,9), 
upon the values of k re i, K c . m , and 9 has been investigated in the case of A=3 and 4. 

In this paper we will proceed further on into this direction by analyzing the contribution of various (ST) channels 
to the one-body momentum distribution of a nucleon Ni belonging to a N1N2 pair in a spin-isospin state (ST). Our 
aim is to understand the quantitative relevance and the momentum dependence of these contributions , in particular 
as far as the deuteron-like state (10) is concerned. In this respect, it should be stressed that in Ref. [Xl|, it has been 
shown that in 3 He and 4 He and for back-to-back nucleons (A^ c m . = 0, the deuteron-like momentum configuration) 
the quantity 

R ^w](krei,K c , m , = 0) = n^™ 0) (fc rcl , K c . m . = Q)/n D (kt e \), (55) 
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i.e. the ratio of the relative momentum distribution of a pn pair in state (ST) — (10) to the deuteron momentum 
distribution, exhibits a constant behavior starting from fc re ; > 1.5 — 2 /to -1 ; this means that at short relative distances, 
the motion of a back-to-back (pn) pair in a nucleus behaves at short distances like in a deuteron. However, a constant 
behavior is not expected to be observed in the ratio of the (ST) one-nucleon momentum distribution to the deuteron 
distribution, since the former, being the integral of the two-body momentum distribution, besides the deuteron-likc 
configuration, includes many other NN configurations. The separation of various (ST) contributions to the one-body 
momentum distribution is an involved task. The problem can be solved by considering the half-diagonal spin-isospin 
dependent two-body density matrix and its integral over r 2 . To this end, it is useful first of all to calculate the number 
of nucleon pairs in a given spin-isospin state. 



A. The number of NN pairs in various spin-isospin states 



The two-body interaction acts differently in states with different spin, isospin and relative orbital momentum L, 
whose values are fixed by the Pauli principle, namely S + T + L = odd. In order to investigate the spin-isospin 
dependence of the momentum distributions, it is useful to start counting the number of pairs N^ 1 ^ in various (ST) 
states in a nucleus with Z protons and N neutrons, with Z + N = A. This quantity is given by Eq. ^ and satisfies 
the sum rule Eq. (fT0|) . The value of Nq£ N2 has been calculated in various papers, e. g. in Refs. [4J,|43] for A < 16, 
in Ref. [23| for A < 4 and in Refs. [461,23] by considering only pairs with L = 0. Here our approach to this topic and 
the results for A = 3, 4, 16, and 40 and L = even and odd will be presented. 

To start with, let us consider a full independent-particle (IP) shell model. In the case of s-shell nuclei the number of 
pairs in (ST) states can readily be obtained. As a matter of fact in A=3 and 4 nuclei the relative orbital momentum 
of all pairs is zero, so that only two (ST) states survive, namely (10) and (01). A pn pair can be either in (10) state, 
with probability 3/4, or in (01) state, with probability 1/4, whereas a pp(nn) pair can only be in (01) state, with 
probability 1. Multiplying these probabilities by the number of pn, pp and nn pairs (NZ, Z(Z — l)/2, N(N — l)/2, 
respectively), the total number of pairs is obtained 

^3(4) = NZ QiV ( t 0) + + ^Zf^-Nffa + ^2^1^ (56) 

from which the total number of pairs in a given (ST) states is obtained, namely 



in 3 He, and 



^( 4 io)=3 Af 01) =3 (58) 



in 4 He (note that the iV( 10 ) state refer to pn pairs only, whereas iV( 01 ) includes pp, nn and pn pairs and it is for this 
reason that no nucleon label appear in NrsT))- In A > 4 nuclei also the states (11) and (00) contribute. In Ref. 
(45j a general approach to calculate, within the IP model, the number of pairs in various (ST) states, based upon 
counting even and odd pairs in spatial configurations corresponding to a given Young tableaux, has been given, and 
explicit formulas can be found there. In our approach the values of N^ T y for the three- and four-nucleon systems 
have been obtained using the wave functions of Ref. [26| and [25[ corresponding to the AV18 and AV8' interaction, 
respectively, whereas for complex nuclei the cluster expansion of Ref. }31j which includes two-, three- and four-body 
cluster contributions has been used to calculate the integral of the diagonal spin-isospin dependent two-body density 
matrix (Eq.©) according to 



N fsT) 



J pfg^ ) 2 (r 1 ,r 2 )drxdr 2 = j ' n^f (k 1 ,k 2 )dk 1 dk 2 = J n^^-(k reU k c . m )dk re idK c . m . (59) 



If IP wave functions are used in Eq. (f5"9")l . the IP values of NA T s have to coincide with the values provided by 
the formulas of Ref. (45|, as indeed it is the case. When the IP model picture is released and a full many-body 
approach with interacting nucleons is considered, odd values of the relative orbital momentum appear also in A=3 
and 4 nuclei so that: (i) the states (00) and (11) are generated in 3 H, 3 He, and 4 He; (ii) the amount of various (ST) 
states in complex nuclei is changed. Thanks to isospin conservation, the number of states (01) is decreased in favor 
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of states (11) and the number of deuteron-like states (10) is also decreased in favor of the state (00). In Ref. |44| 
^(st) 2 nas been calculated for 3 He, 4 He, 6 Li, 7 Li, and 16 O using variational Green's Functions Monte Carlo wave 



functions and various Argonne interactions; in Ref. 23[ Nfy T -, has been obtained for nuclei 3 He, 3 H and 4 He using 
wave functions resulting from the correlated Gaussian basis approach [32| and the V8' interaction, finally, in Ref. 
(46l |47| the number of pairs in L=0 state has been evaluated through the Periodic Table using phenomenological 
correlated wave functions. In the present paper we have calculated N^ T ~, for A — 3, 4, 16, 40 using wave functions 
obtained within the hyperspherical harmonic variational method [26| and the AF18 interaction, for A = 3, the ATMS 
method of Ref. [24l l25j and the AV8' interaction, for A = 4, the linked-cluster expansion of Ref. [3l[ and the AV8', 
for A = 16 and A — 40. The results of our calculations, which are presented in Table IIII1 clearly show that: (i) 
there is satisfactory general agreement between our results and the ones of Ref. [23j, |44j ; (ii) as previously found 
in those papers, when the IP model picture is released and NN correlations are taken into account, the value of 
-^(10) is practically unchanged, whereas the number of pairs in the (01) state is decreased in favor of the state (11). 
The reason for that was nicely explained in Ref. [23|, |44| : it is due to some kind of many-body effects induced by 
tensor correlations between particles " 2" and " 3" , generating a spin flip of particle " 2" , and giving rise to the state 
(11) between particles "2" and "1". These effects are automatically included in our calculations, since "exact" wave 
functions are used in case of few-nucleon systems and a cluster expansion embodying many-body clusters is adopted 
in our approach for complex nuclei. 

B. The spin-isospin contributions to the momentum distribution 

We apply here Eq. (|26p (with k\ = k), obtaining for the proton momentum distributions in 3 He 

n P 3 (k) = ln p T=0 (k) + ^n p T=1 (k)= (60) 
= nf°\k) + nT°\k) + nT 1 \k) + nf 1 \k), (61) 
since there is only one pp and one pn pair containing proton " 1" , whereas the neutron distribution is given by 

nm = \n n T=0 (k) + -n n T=1 {k) = (62) 

= ^ l(10) (fc) + ^ (00) (fc)+^ (01) (fc)+^ (11) (fc), (63) 

since there are two pn pairs containing neutron "1" and no pp pairs. The momentum distributions of 4 He, le O and 
40 Ca are given, respectively, by 

u 4 (fe) = \n [ ^l{k)+ l -n T=1 {k)= (64) 

= n^W+nf^W+nf^W+n^ik), (65) 

n 16 (fei) = ^:l(k) + ^n T=1 (k) = (66) 
= ng°\k) + n%°\k) + n[f\k) + nT\k), (67) 

n 40 (fc) = l n «(fe) + gn T=1 (fe)= (68) 

= ng 0) (k) + 4° 0) (fc) + (k) + n$ ll) (k), (69) 
where Eqs. (|23|) and ([22| have been used, 

4 P :l(k)=4 P =i(k)^T=i(k), (70) 

and 

J n A {k)dk = J n%l (k)dk = J n T =i{k) dk = I. (71) 
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The results of calculations of the spin-isospin contributions to the momentum distribution of 3 He, 4 He, 16 and 40 Ca, 
are presented in Figs. [TTIITS1 the following remarks are in order: (i) the contribution from the (00) state is negligible 
both in few-nucleon systems and complex nuclei; (ii) the (11) state in 3 He and 4 He is small, both at low and large 
values of k, but it plays a relevant role in the region 1.5 ;$ k < 3/m _1 ; (iii) in the proton distribution of 3 He (Fig. 
Hip the (01) contribution is important everywhere except in the region 1.5 < k < whereas in the neutron 

distributions (Fig. [T2]), thanks to the different weight of the (01) state (1/4 instead of 5/8; cf. Eqs. (|B(J)) and ([52])). 
the contribution from this state is much smaller; (iv) in complex nuclei the (11) state (odd relative orbital momenta) 
plays a dominant role, both in the independent particle model and in the many-body approach (cf. Table IIIII and 
Figs. 1141 and [TS"]) . Thus, in summary, we found that all spin-isospin components, except the (00) one, contribute to 
the high momentum content of the momentum distributions and only in the case of the neutron distribution in the 
non iso-scalar nucleus 3 He, the deuteron-like state (10) is the dominant contribution. 

In order to provide further evidence of the universality of SRC, we show in Fig. [16] the "elementary" quantities 
n T™o(^i) an d «T=i(fei) for different nuclei, and it can be seen that, starting from k = \ki\ ~ 2fm~ 1 , they follow 
the same pattern. 



V. THE MOMENTUM DISTRIBUTIONS OF NUCLEI VS THE DEUTERON MOMENTUM 

DISTRIBUTION 

As it clearly appears in Fig. [51 at k > 1.5 — 2/m _1 the momentum distribution of nuclei exhibits a trend similar 
to the one of the deuteron. However a quantitative analysis of the ratio 



nA(k) 
n D (k) 



RA, D {k) = ^ (72) 



is in order, since nAik) is usually interpreted as the scaled deuteron momentum distribution, i.e. Ra/dQz) = 
nA(k)/riD(k) ~ const. Such an interpretation originated long ago either from the use of pioneering theoretical 
many-body calculations [T7i - f20| . or by assuming it as an input for the calculations of n^(fci) at k > kp plj when 
variational Monte Carlo calculations were difficult to perform at high values of the momentum, or by obtaining the 
momentum distributions from an average value of the pn and pp spectral functions [l2| . Having nowadays at disposal 
advanced many-body calculations of the momentum distributions performed with realistic models of the two-nuclcon 
interactions, a quantitative analysis of Eq. ([72]) is timing. To this end, we show in Fig. [17] the ratio i?A/r>(^) cal- 
culated with realistic many-body wave functions. It clearly appears that starting from k > 2/to~ 1 , the ratio is not 
a constant but appreciably increases with k. Let us discuss the origin of such an increase. A first possible origin 
should be sought in the different role played by pn and pp correlations. As a matter of fact, the proton and neutron 
momentum distributions in 3 He shown in Fig. 1181 exhibit a different rate of increase, which can qualitatively be 
understood in terms of SRC as follows: in 3 He the proton momentum distribution is affected by SRC acting in one 
pn and one pp pairs, in the former pair the deuteron-like state (10) is three times larger than the (01) state, whereas 
in the latter pair the deuteron-like state is totally missing; on the contrary, the neutron distribution is affected by 
SRC acting in two proton-neutron pairs, with a pronounced dominance of the deuteron-like state (10); therefore, one 
expects that around k ~ 2/m~ 1 , where, np SRC dominate over pp SRC [TUHl]), n ^l n D — 2 and n^/un — 1, which 
indeed seems to be the case. However, other effects of different origin can contribute to the deviation of the ratio 
nA(k)/no(k) from a constant. These are due to the center-of-mass motion of a pn pair in a nucleus, to the different 
role played by the states (01) and (11) in different nuclei, and, particularly, to the fact that being the one-nucleon 
momentum distribution the integral of the two-body distribution over fe, rtyi(ki) = J n^(ki, k2)e?fc2, it may contain 
configurations different from the deuteron one (back-to-back nucleons) . In order to better investigate these possibility, 
let us consider the spin-isospin ratio 



n D (k) 



«J52(*) = ^# (73) 



which is shown in Figs. IT5][2"51 (note that, as mentioned in the caption of Fig. [121 the quantity n^f T ' includes the 
proper coefficients which multiply the "elementary" quantities n^ 1 * ). It can be seen that the behavior of the proton 
and neutron ratios for 3 He clearly shows that in the region 1.5 < k < 3 fm _1 the former is governed by the (01) state 
in the pp and pn pairs; on the contrary, the neutron ratio is fully dominated by the deuteron-like (10) state in the two 
pn pairs. The most interesting ratio is R^^ D (k) — (fc)/ri£>(fe), since it provides information on the behavior of 

the deuteron-like pairs in nuclei; it can be seen that in the region of SRC (k > 2fm~ 1 ) R ( ™ ] D increases with increasing 
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value of k, with a different rate of increase for different nuclei: it is about 30% in the neutron momentum distribution 
of 3 He, and of the order of 100% in other nuclei. As already pointed out, the increase of the ratio R^J^{k) with k 
could also be due to the cm. motion of the pair in the nucleus. To take this into account is no easy task. As a matter 
of fact, consider the simple case when the (10) two-body momentum distribution factorizes in the following form [l2j 

n p ^ 0) (ki,k 2 ) = nfo (k re l, K c . m ) = n P Q {krei, K c . m ., 0) ~ n D (k re i)n^ \ m {K c . m .) , (74) 

where n c .m,{K c . m ,), calculated from a many-body approach in Ref. [ill ], can be approximated by a OS* wave function. 
In Refs. [llj and [48j . Eq. ([7S|) has indeed been shown to hold, but only in a restricted region of k re i and K c , m ,, 
namely 

K c . m < 1.0 - 2.0 fm- 1 k rel > k; el = f A (K c . m .) (75) 

where the function f A depends upon K c m and A in such a way that the value of k~ el increases with increasing values 
of Kcm- Thanks to momentum conservation k 2 — — (fei + K c . m ), one can write 

n (To)(^ 1 ) — / n D(\ki — \)n c (\K c . m .\)dK c . m . (76) 

which shows that only in the case of a pn pair at rest, i.e. n c . m .(K cm.) = 5(K c . m .), one has ^fo(^i) — n D{k\) 
R a /jj (fci) — const. The convolution of the deuteron momentum distributions with the cm. motion leads to an increase 

of r^Q l (fci), whose magnitude and rate of increase depend upon the detailed forms of n\^{k re i) and n c-m .(K c . m ,); 
moreover, since, as already stressed, the one-body momentum distribution is the integral of the two-body momentum 
distribution, configurations different from the factorized one (Eq. (17410 can contribute to the integral (Eq. (|76p). 

A. On the short-range deuteron-like configurations in nuclei 

A particular useful quantity to understand SRC in nuclei is the one that is obtained by integrating the two-nucleon 
momentum distribution of the state (10) in a narrow range of the cm. momentum (K c .m. $ 1 — 1.5 /to -1 ), when the 
cm. and relative motions are decoupled, and Eq. (j?4")) is satisfied [llj, namely 



n P p /A (k re i) = J n p " 0) (k re i,K c . m ., 6) d K c .m. ^ n D (k rei ) 4tt J n c m {K c . m )K c m dK c . m . (77) 

In Ref. the 2N SRC probability in the deuteron has been defined as the integral of the deuteron momentum 
distribution in the range k re i > 1.5 fm -1 (cf. Table HV)) . therefore we can consider as the analog in a nucleus the 
quantity 



^_D/A=47r/ n p ™ /A (k re i)k rel dk rel (78) 
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where n p ^ A is given by Eq. (f77j) . We can also define the total number of quasi-deuteron short-range correlated pairs 
as follows 

N d ,a = Nf 10) V D/A (79) 

where the number of N^ Q ^ pairs is listed in Table IIIII The calculated values of Vjj /a an d N D /a are given in Tables 
IIVI and |V] The decreasing behavior of Njj /a with A can easily be understood as due to the increasing importance 
of higher momentum components in the cm. momentum distribution of the pair when A increases, with a resulting 
flatter distributions in heavier nuclei (cf. Fig. [2~4l ; as a result only a smaller part of the distribution is included in 
the integral over i4T c m .. We have also considered the quantity a 2 [J], the per-nucleon probability of deuteron-like 
configurations in A with respect to the probability of SRC in the deuteron (~ 0.04). Our values for A < 40 are less 
than the values extracted from the A(e,e')X experiments however such a comparison is perhaps a premature 

one, since, from one side, other deuteron-like configurations which occur outside the factorization region should be 
considered in the theoretical calculation, and, from the other side, a careful investigation of the effects of FSI effects 
on the extraction of a 2 from the inclusive A(e, e')X cross section ratio might be in order. We should also mention, in 
this respect, that the values of a 2 were also recently calculated in ref. [46[ within an approach in which only correlated 
pairs in L — states are considered (cf. Table IIIII) , obtaining results that coincide with the ones obtained in the 
present paper for A = 3, 4, and are a bit lower for A > 4. 
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VI. SUMMARY AND CONCLUSIONS 



Recently, the universal character of SRC in few-nucleon systems ( 2 H, 3 H, 3 He, and 4 He) has been demonstrated 
by calculating the dependence of two-body momentum distributions upon the relative momentum |k re /| = k re i of the 
correlated pair (23|, as well as upon the center-of-mass momentum |K c m . | = K c . m . and the angle between K cm and k re ; 



In the present paper we have addressed the problem of the effects of SRC, and their spin-isospin components, 



on the one-nucleon momentum distributions riA(k) of few-nucleon systems and complex nuclei. The momentum 
distribution is per se a relevant quantity in nuclear theory, moreover, it plays a relevant role in the interpretation 
of various experimental data, in particular in inclusive experiments of lepton scattering off nuclei at medium and 
high energies. Using the proper diagonal and non diagonal one- and two-body spin- and isospin- dependent density 
matrices, we derived in Section [TT] the expression of the momentum distributions of a nucleon belonging to a NN pair 
in a state with total spin S and isospin T. Before calculating the spin-isospin structure of the momentum distributions, 
in Section inil we presented some general concepts concerning nucleon momentum distributions and a clear cut way to 
separate them in a mean-field and SRC contributions, and after that we have analyzed the results of the most recent 
calculations of the momentum distributions for nuclei with A=2, 3, 4, 16 and 40, performed within realistic many-body 
approaches and modern NN interactions. The aim was to ascertain whether some general features of the momentum 
distributions could be established, which are independent, to a large extent, of the NN interaction and many-body 
approach. The results of our analysis show indeed that, even if quantitative differences are provided by different 
interactions and many-body approaches, some general and universal features of the momentum distributions, can be 
single out, namely: (i) at k < 1 — 1.5 fm , the mean-field approach dominates the distributions, with a resulting 
sizeable A dependence; (ii) at larger values of k, of the order of 2/to~ 1 , due to the effects of SRC, the momentum 
distributions abruptly change their slope, exhibiting a k dependence which is very similar in different nuclei, hereby 
demonstrating the universal character of SRC; (iii) the correlated part of the momentum distribution is by orders 
of magnitude larger than the predictions of any mean-field approach, so that experiments providing even rough 
information on high momentum components would be able to rule out mean-field predictions. Similar conclusions, 
reached in the past by phenomenological calculations (see e.g. Refs. [12] and [HI), are quantitatively confirmed by 
the present systematic analysis. After checking that the calculation of the high momentum part of riA{k) is well 
under control, we turned in Section ITVl to the calculation of the spin-isospin structure of the momentum distributions. 
First of all we calculated the number of NN pairs in various spin-isospin states in different nuclei, both within the 
independent particle models and in many-body approaches embodying SRC, finding agreement with calculations 
performed within different approaches and confirming that SRC have very small effects on the number of isosinglet 
pairs in state (10), unlike what happens with isotriplet pairs in state (01), whose number is decreased in favor of the 
pairs in (11) state. We have calculated the contribution of the states (ST)=(10), (00), (01) and (11) to the momentum 
distributions, finding that all of them, except the state (00), have comparable effects in a wide range of momentum. 
The contribution of the isosinglet state T=0 is almost entirely exhausted by the (10) state, whereas both states (01) 
and (11) contribute to the isotriplet state (11). We found that at momentum values k > 2/m -1 , the contribution of 
both isosinglet and isotriplet states follow the same pattern, which represents further evidence of the universality of 
SRC. A systematic and quantitative comparison between riA(k), and its spin-isospin components, n^ T (fc), with the 
deuteron momentum distribution, n£>(k), has been presented in section HVl by analyzing the ratios n^ T (fc)/n£)(fc). 
We found that in the region of SRC, k > 2/m -1 , this ratio does not stay constant but increases with increasing 
k, and interpreted such a behavior as due to the presence in the momentum distribution of two-nucleon momentum 
configurations arising from the cm. motion of a pair and differing from the back-to-back nucleons configuration. Our 
spin-isospin dependent approach allowed us to calculate also: (i) the relative momentum distribution of a proton- 
neutron pair moving with small cm. momentum and its integral in the range 1.5 < k < oo, a quantity which is 
assumed to represents the probability of two nucleon SRC in a nucleus, finding similar values (~ 0.04) in a wide range 
of A, 2 < A < 40; (ii) the total number of SRC pairs in (10) state, interpreting its A-dependence in terms of the 
A-dependence of the cm. momentum distribution; (iii) the per-nucleon probability of two-nucleon SRC in nuclei, a 
quantity which is under active experimental investigation. To sum up, we would like to stress that the theoretical and 
experimental investigations of SRC in nuclei is progressing very rapidly (see e.g. [U, E3) and we do expect that the 
results presented in this paper will contribute to clarify several still open problems in this important field of hadronic 
physics. 
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MEAN FIELD AND CORRELATION PROBABILITIES 


Nucleus 


Potential 


Pgr 




3 He [26] 


AV18 [39] 


0.677 


0.323 


4 He [21, 24] 


RSC [35] AV8'[37] 


0.8 


0.2 


Nucleus 


Potential 


Po 


Pi 


16 [22] 


V8' [37] 


0.8 


0.2 


40 Ca [11] 


V8' [37] 


0.8 


0.2 



TABLE I: The proton probabilities V p gr(0) = / dki n p r(0) (kx) and P p ex{1) = / dk x (fei) 





2 H 


3 Hc(n) 


3 He(p) 


4 He 


16Q 


40 Ca 


fcf [fm" 1 ] 


P 


Pi 


V 


Pi 


Vo 


Vi 


V 


Vi 


Po 


Pi 


0.00 


1.000 


0.999 


0.677 


0.323 


0.84621 


0.15285 


0.79999 


0.20016 


0.80 


0.19321 


0.50 


0.3078 


0.568 


0.277 


0.201 


0.53643 


0.14032 


0.66972 


0.19635 


0.69997 


0.18301 


1.00 


0.081 


0.163 


0.038 


0.0723 


0.10479 


0.1045 


0.17588 


0.14794 


0.24706 


0.13771 


1.50 


0.0366 


0.067 


0.0049 


0.036 


0.0079 


0.0791 


0.00792 


0.09417 


0.01022 


0.10143 


2.00 


0.0221 


0.041 


0.0015 


0.024 


6.9512 10" 4 


0.06156 


5.9 10" 5 


0.06344 


3.28 10" 4 


0.07124 



TABLE II: The values of the partial probability, Eq. (H5]l . for 3 He, 4 He, and 16 and 40 Ca, calculated for different values of 
the momentum fc^ with = oo, . 
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1.50 
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Me 


oKt (.Present work) 


1.488 


1.360 


0.013 
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bKC 1441 


i k a 


i 9 a 


a ai 
U.U1 


A 1 A 
U.14 




CDP Tool 
OXlO |ZO] 


1.489 


1.361 


0.011 


0.139 




IPM 


3 


3 








Yr\slY\\}S StatCSJ |4DJ 


3 


3 






-tie 


oxio ^.r resent worxj 


o on 

z.yy 




A A1 
U.U1 


A A Q 




.3 1 \ 1 | ^L^iJ 


3.02 


2.5 


0.01 


0.47 




nnri Tool 

bKC |23| 


o aaq 


2.01 2 


A AO 

0.08 


A /tOO 

0.428 




IP IVi 


30 


30 


6 


54 




IPM(0s states) [46] 


20 


18 






16 


SRC (Present work) 


29.8 


27.5 


6.075 


56.7 




SRC [44] 


30.05 


28.4 


6.05 


55.5 




IPM 


165 


165 


45 


405 


40 Ca 


IPM(0s states) [46] 


90 


20 








SRC (Present work) 


165.18 


159.39 


45.10 


410.34 



TABLE III: The number of pairs N^j,^, Eq. ([9]), (only pairs in relative L = motion in Ref. [4(|) in various spin-isospin states 
in the independent particle model (IPM) and taking into account SRC within many-body theories with realistic interactions. 



2N SRC PROBABILITY 


N N Interaction 


2 H 


3 He(n) 


3 He(p) 


4 He 


i6 


40 Ca 


RSC 


0.04 






0.09 


0.12 




AV14 


0.036 






0.11 


0.14 




AV8' 


0.036 






0.09 


0.10 


0.10 


AV18 


0.037 


0.067 


0.041 








CS 


0.033 


0.079 


0.046 


0.09 


0.10 


0.14 



TABLE IV: The value of the two-nucleon short-range correlation (2N SRC) probability in the deuteron, 4n no(fc) k 2 dk, 
and in complex nuclei (Eq. (|78|ll obtained with momentum distribution resulting from many-body calculations performed with 
different NN interactions. The result (CS) of the phenomenological model of Ref. [T3] is also shown. 





2 H 


3 He 


4 He 


16Q 


40 Ca 




Vd/a 


N D/A 




Vd/a 


N D/A 




Vd/a 


N D /A 


a 2 


V D /A 


N D /A 


0(2 


Vd/a 


N D /A 




1.5 


0.04 


0.04 


1 


0.04 


0.06 


1 


0.04 


0.12 


1.5 


0.031 


0.93 


2.9 


0.030 


4.9 


6.1 



TABLE V: The values of V D /a (Eq. (J78J)) and N D /A calculated in correspondence of K£ m . = 1.5 fm 1 . The quantity 

ci2 — [(2/A)][N D / A /N D / D ] is the per-nucleon probability of 2N SRC in A with respect to the deuteron. 
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k [fin 1 ] k [fin" 1 ] 

FIG. 1: (Color online) Deuteron momentum distributions in logarithmic (a) and linear (b) scales corresponding to various NN 
interactions: RSC [HI], Paris |3fJ, AV8' H3], AV14 [H| and AV18 [H]. Unless otherwise stated, here, and in the other Figures, 
the normalization is 47r J k 2 dkn ] ^ 1 (k) = 1. In this and the following Figures |ki| = k. 
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k [fin 1 ] k [fm" 1 ] 

FIG. 2: (Color online) The proton and neutron momentum distributions of 3 He in logarithmic (a) and linear (b) scales. Three- 
nucleon wave functions from Ref. [2(|. The full curve represents the deuteron momentum distribution. Both 3 He and deuteron 
wave functions correspond to the AV18 interaction [3^ ]. 
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k [fin -1 ] 
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2.5 3.0 3.5 

k [fin ] 



4.0 



FIG. 4: (Color online) The momentum distribution of 16 in logarithmic (a) and linear (b) scales, corresponding to different 
wave functions and NN interactions: Dash: Ref. [Tt]]; Dots: Ref. [%|; dot-dash: Ref. [3l]]. The parametrization of Ref. [12] is 
also shown by the short-dashed curve (CS). The full curve represents the deuteron momentum distribution corresponding to 
the AV18 interaction. 
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FIG. 5: (Color online) The momentum distribution of 40 Ca in logarithmic (a) and linear (b) scales, corresponding to the AV8' NN 
interaction and calculated within two different many-body frameworks. Dash: cluster expansion up to FHNC/SOC order [2^]; 
Dots: cluster expansion at second order 31]. The full curve represents the deuteron momentum distribution corresponding to 
the AVI8 interaction. 
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FIG. 6: (Color online) The proton momentum distribution of nuclei considered in this work in logarithmic (a) and linear (b) 
scales, calculated within different many-body approaches with equivalent NN interactions, namely the AV18 one, in the case 
of 2 H and 3 He, and the AV8' one, in the case of 4 He, 16 0, and 40 Ca. 
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20 




FIG. 9: (Color online) The same as in Fig. [7] but for ie O. Wave functions from Ref. [U, AV8' interaction. The values of 
Vo = 47T J k 2 dk no (k) and Vi = 4tt / fc 2 dfc ni are listed in Table HI 
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k [fin 1 ] 

FIG. 11: (Color online) The various spin-isospin contributions, (ST), to the proton momentum distribution of 3 He. Wave 
function from Ref. [26j |. AV18 interaction. The continous line without symbols is the sum of the four contributions (cf. Eq. 

(ED). 
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FIG. 12: (Color online) The same as in Fig. [TT]but for the neutron distribution (cf. Eq. (|63(l \ 
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FIG. 13: (Color online) The various spin-isospin contributions to the proton momentum distribution of 4 He (cf. Eq. I)65| l), 
Wave function from Ref. [25j], AV8' interaction. 
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FIG. 17: (Color online) The ratio of the proton momentum distribution of nucleus A shown in the previous Figures to the 
deuteron momentum distributions. 
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n,(k)/n n (k) 



n-"'(k)/n D (k) - * - "(k)/n D (k) 




1 2.3 

k [fm ] 



FIG. 19: (Color online) The various spin-isospin contributions to the ratio of the proton momentum distributions of 3 He (Eq. 
(jGlJ) to the deuteron momentum distributions. Wave function from Ref. [261 ] . AV18 interaction. 
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FIG. 21: (Color online) The various spin-isospin contributions to the ratio of the proton momentum distributions of 4 He to the 
deuteron momentum distributions. Wave function from Ref. [2(|, AV8' interaction. 
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FIG. 22: (Color online) The various spin-isospin contributions to the ratio of the momentum distribution of O to the deuteron 
momentum distributions. Wave function from Ref. [31|. AV8' interaction. 
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FIG. 23: (Color online) The various spin-isospin contributions to the ratio of the momentum distribution of Ca to the 
deuteron momentum distributions. Wave function from Ref. [3l[, AV8' interaction. 
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